Quantum tunneling of magnetization in single molecular magnets coupled to 

ferromagnetic reservoirs 
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The role of spin polarized reservoirs in quantum tunneling of magnetization and relaxation pro- 
cesses in a single molecular magnet (SMM) is investigated theoretically. The SMM is exchange- 
coupled to the reservoirs and also subjected to a magnetic field varying in time, which enables 
the quantum tunneling of magnetization. The spin relaxation times are calculated from the Fermi 
golden rule. The exchange interaction of SMM and electrons in the leads is shown to affect the spin 
reversal due to quantum tunneling of magnetization. It is shown that the switching is associated 
with transfer of a certain charge between the leads. 
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Introduction Single molecular magnets (SMMs) form 
a class of systems whose permanent magnetic moments 
stem from their molecular structure [H, Generally, 
SMMs are characterized by a large ground state spin 
number S and a relatively large uniaxial-type magnetic 
anisotropy. As a result, an energy barrier appears for 
switching SMM's spin between the two stable spin states 
| ± S). At higher temperatures SMMs behave like param- 
agnetic or superparamagnetic particles with a large mag- 
netic moment. When temperature is lowered, the ther- 
mal energy is not sufficient to reverse spin orientation of 
the molecule. 

It has been recently predicted that the molecule's spin 
can be reversed by a spin current 0, 0] • This method al- 
lows magnetic switching without any external magnetic 
field and is of great interest from the point of view of 
future applications in spintronics and information tech- 
nology [5|, |6| . Another way to switch spin of the molecule 
relies on the phenomenon of quantum tunneling of mag- 
netization (QTM) 0], which occurs in a magnetic field 
varying (linearly) in time. The quantum tunneling takes 
then place between the energetically matched levels on 
the opposite sides of the barrier, and leads to successive 
steps observed in hysteresis loops [||. 

An important problem associated with the phe- 
nomenon of QTM is the spin relaxation. Such a relax- 
ation can take place for instance via spin-orbit or spin- 
spin couplings. In this paper we consider spin switching 
via QTM, associated with the spin relaxation due to cou- 
pling of the molecule with two reservoirs of spin polarized 
electrons. More specifically, the system under consider- 
ation consists of a SMM located in a junction, and sep- 
arated from two ferromagnetic leads by tunnel barriers, 
see fig. [lja). For simplicity, we take into account only 
collinear (parallel and antiparallel) magnetic configura- 
tions of the leads. Moreover, the magnetic easy axis of 
the molecule is parallel to the magnetic moments of the 
leads. To allow the QTM, an external time-dependent 
magnetic field is applied to the molecule along its easy 
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FIG. 1: (color online) (a) Schema of the system under consid- 
eration for parallel (black solid arrows) and antiparallel (grey 
solid arrows) magnetic configurations, (b) Energy levels cor- 
responding to different spin states of the SMM. The grey dots 
represent the occupancy of the molecular spin states when the 
third resonant field is reached. 



axis, sec fig. H^b). We assume that the field has no 
influence on the magnetic configuration of the system, 
which can be controlled by different means (exchange 
anisotropy, for instance). In addition, we assume a signif- 
icantly weaker magnetic field perpendicular to the mag- 
netic easy axis. 

Theoretical description In a recent paper Q we con- 
sidered the situation when a bias voltage was applied be- 
tween the two external ferromagnetic leads, and the spin 
reversal was due to the associated spin current. Here, 
we consider a different situation, i.e. the system is un- 
biased but instead an external magnetic field is applied 
to switch the molecule owing to the QTM phenomenon. 
The switching is accompanied by a current impulse be- 
tween the leads, which is a consequence of the spin relax- 
ation processes. This, in turn, is associated with transfer 
of a certain charge between the leads. 

The simplest model Hamiltonian describing the SMM 
takes the form HflSuUl 



MM 



7~Lo + Ti-QTM + Hi. 



(1) 
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The first term includes the uniaxial anisotropy and the 
Zeeman energy of the molecule in a field along the 
anisotropy axis (the field varies linearly in time), 



Ho = -DS 2 Z + g[i B S z H z 



(2) 



where S z is the z-component of the spin operator, D is 
the uniaxial magnetic anisotropy constant, g is the Lande 
factor, and fis stands for the Bohr magneton. 

The next term, Hqtm, of the Hamiltonian is involved 
in the phenomenon of QTM, and includes the transverse 
anisotropy and the Zeeman energy of the molecule in a 
field perpendicular to the anisotropy axis, 



Hqtm = E(S X - Sy 
+ 9Hb{S x H x 



f C(S% - 

SyHy), 



st 



(3) 



where S x and S y are the transverse components of the 
spin operator, S± = S x ± iS y , whereas E and C are the 
transverse magnetic anisotropy constants. The perpen- 
dicular magnetic field is much weaker than the longitudi- 
nal component H z , and is assumed to grow in time with 
the same time scale as H z , i.e. H x = 0.1H Z (H y = for 
simplicity) . 

The last term of the Hamiltonian, Hi, describes ex- 
change interaction between electrons of both leads and 
the SMM HHHm, 



Hi 
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J, 



H.c. (4) 



Here, S is the SMM's spin operator, cr = (a x , a y , <j z ) 
is the Pauli spin operator for conduction electrons, a? 
( a kL) are ^ ne annihilation (creation) operators of elec- 
trons in the left (q = L) or right (q = R) electrodes, k (k') 
denotes a wave vector, and a ((3) is the electron spin in- 
dex. In eq. (j4|) J q q i is the relevant exchange parameter, 
assumed to be independent of energy and polarization of 
the leads. In a general case Jl,l ^ Jr,r. ^= Jl.r = Jr,l- 
In the following, however, we assume the symmetrical 
situation, where Jl,l — Jr.r — Jl,r — Jr.l = J- Ow- 
ing to the proper normalization, J is also independent of 
the electodes' size. Finally, in eq. (@| N q (q = L, R) de- 
notes the number of elementary cells in the q-th electrode. 
Since we consider only the unbiased situation, we have 
omitted the direct tunneling between the leads (i.e., tun- 
neling without exchange interaction with the molecule). 

First, we present some general equations describing 
spin relaxation of the molecule, and start from the tran- 
„>(<) 
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7m + 7„ 



LL>(<) 
7m 



sition rates 7„ 

RR > ( <C ) 

7m , at which the molecule's spin changes from 
m) to the upper (lower) state \m ± 1) due to cou- 
pling to the leads. Here, 7^^^ is the transition rate 
of electrons tunneling from the left to right electrode, 
whereas 7m i> ^ < ^ corresponds to electrons tunneling the 
other way. Similar meaning have the rates 7 m L> ^ < ^ and 



7m R> ^ < ^- They refer, however, to electrons which when 
interacting with the molecule change its magnetic state, 
but remain in the same electrode. All these transition 
rates can be calculated from the Fermi golden rule, and 
for instance 



7 ; 



^f >(<) = E E ^m ± l f(4a) [1 - f(4 P )] . (5) 



kct k'/3 



where /(e) is the Fermi-Dirac distribution, and e£ Q (q — 
L,R) is the energy of electron states in the leads. The 
probability of electron transitions from the initial state 
i <-> {ikam} to the final one j «-> {i?k'/3m ± 1} is given 
by W] = {2Tr/h)\{j\H I \i)\ 2 6(E j - E t ), where E t = + 
E m , and E rn is the energy of the unperturbed molecular 
spin state \m). Energy Ej of the final state is given by 
a similar expression. Finally, similar formulae also hold 
for the transition rates 7 m L> ^ < \ 7m L> ^ and 7™ 



The final expressions for the transition rates 7, 
given by the formulae 
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where D q G is the density of states (DOS) at the Fermi 
level in the g-th electrode for spin a, A±(m) is defined 
as A ± (m) = S(S + 1) - m(m ± 1), and C(e) = e[l - 
exp(-e/3)] _1 with /3 _1 = k B T. 

Since our objective is to calculate the average value 
of the z component of SMM's spin, (S z ) — J2 m iTiP m , 
we have to determine the probabilities P m of finding the 
SMM in the spin state \m). To start with, we assume 
the initial spin state of the SMM to be | — S). Taking 
into account negative sign of the gyromagnetic factor, we 
write H z = —H. Thus, in the initial state the magnetic 
field H is negative. The field H grows then linearly in 
time and the SMM's spin undergoes transitions (due to 
QTM) at relevant resonant fields from the state | — S) to 
states IS" — M) (consecutively for M = 0, M = 1, etc) on 
the opposite side of the energy barrier. The probabilities 
of all other states (for each M), i.e. of states | — S + 
l),...,\S — M — l), are then equal to zero. If the SMM's 
spin tunnels to the state other than \S), the interaction of 
the molecule with the electron reservoirs leads to further 
relaxation of the spin to the final state \S), fig. HJb). 

There are two time scales set respectively by the speed 
c at which the magnetic field is increased, c = dH/dt, 
and the relaxation processes due to interaction with the 
electrodes, with the latter scale being much shorter. This 
fact can be used to simplify the problem, and the proba- 
bilities P m can be then found from the set of master equa- 
tions, separately for each field range between the succes- 
sive resonant fields at which the QTM occurs. More- 
over, when TcD, the relevant relaxation processes are 



3 



those which lower energy of the molecule (and increase 
the quantum number m in the case under consideration) . 
For the M-th range the equations take then the form 



(cPs-M 
C Prn 
CPs 



~ 7m + 7m- 
71-1^5-1, 



1 Prr 



(7) 



for m e (S-M+l, 5-1), and P defined as P = dP/dH. 
One should bear in mind that the transition rates 7 m 
depend on magnetic field (see eq. ©). Additionally, we 
note that 7 m « for m > 1, which justifies the absence of 
terms corresponding to the transitions to lower molecular 
spin states in eqs. (0. 

The boundary conditions for the probabilities in 
eqs. J7]) (at the resonant field H^) are: 



Ps-M\fl 



(M)\ _ 



— Ps-M, 

pM-ll H {M) 



(8) 



for m € (5- M + 1,5). In eqs. ® P 7 l(H (n +^) de- 
notes the probability P m calculated for the field range 
{#(«), #(«+*)) and taken at the field H i ~ n+1 \ More- 
over, Ps-n is the probability of the QTM to occur 
between the states | — S) and |5 — n). This prob- 
ability can be obtained analytically with the repeti- 
tive use of the two-level Landau-Zener model [3, [HI], 
Ps-n = (l - exp[-27rA„]) H™~ exp[-27rAj], where A„ = 
A 2 l /[Ah(2S — n)gfj,Bc\ and A n is the splitting of the two 
states due to the TLqtm term in Hamiltonian. 

Numerical results Numerical calculations have been 
obtained for an octanuclcar iron(III) oxo-hydroxo clus- 
ter of the formula [Fes02(OH)i 2 (tacn) 6 ] 8+ (shortly Fes), 
whose total spin number is S = 10. The anizotropy con- 
stants D = 0.292 K, E = 0.046 K, C = -2.9 x 10" 5 K as 
well as the splittings A n are adopted from refs 
We assume that J w 100 mcV. Furthermore, 



11, 1 



for both 

leads we assume the elementary cells are occupied by 
2 atoms contributing 2 electrons each. The density of 



free electrons is assumed to be n ~ 10 



29, 



The elec- 



trodes are characterized by the polarization parameter 
Pi = (D q + - D q _)/(D\ + D q _), where D q +{ _ } denotes 
the DOS of majority (minority) electrons in the g-th elec- 
trode. The temperature of the system is assumed to be 
T = 0.01 K, which is below the blocking temperature 
T B = 0.36 K of Fe 8 . 

Let us consider first the average value of the z com- 
ponent of SMM's spin, (5 Z ), in an external magnetic 
field increasing linearly in time, fig. [2j The reversal of 
the SMM's spin is due to the QTM, whose signature is 
the characteristic staircase pattern. The steps occur at 
the resonant fields where the QTM is allowed, and their 
heights are determined by the field sweeping speed c and 
parameters of the TLqtm term. For small values of c 
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FIG. 2: (color online) The average value of the SMM's spin, 
{Sz}, as a function of external magnetic field H for various 
field sweeping speeds c. The magnetic moments of the elec- 
trodes are in the parallel configuration and P L = P R = 0.5. 
Solid lines represent the case of the SMM's spin coupled to 
electrons, whereas dashed ones describe the system without 
the exchange interaction between the molecule and electrons 
(and in the absence of other spin relaxation processes). 



full reversal occurs already at the third resonance field 
(corresponding to n = 2). The exchange coupling with 
reservoirs is essential to observe the full reversal of the 
spin - after QTM the spin relaxes to the state \S). 

Typical values of the transition times l/7 m (for the 
parameters assumed) are of the order of 10" 13 — 10" 12 
s. The duration of the SMM's spin reversal is thus set 
by the rate c at which the field H is augmented, and it 
varies between 1 s, fig. [2jd), and 10 3 s, fig[2fa). As a 
consequence, the steps in figs. [2] are very sharp. 

Generally, the values of 7^ depend on the magnetic 
configuration of the system. Assume the same total den- 
sity of states in both electrodes, but generally differ- 
ent spin polarizations. In the parallel (P) configuration 
7^J > oc 1 — (P L + P R ) 2 /4, whereas in the antiparallel 
(AP) one 7^ p> « 1 - (P L - P R ) 2 /A. If additionally 
we assume now P L = P R = P, the above formulae re- 
duce to 7^ > cx 1 — P 2 for the parallel configuration, 
while in the antiparallel configuration 7^ > is indepen- 
dent on the polarization P. The difference between these 
two configurations can be clearly seen for P = 1, which 
corresponds to the situation with both electrodes being 
perfect halfmetallic ferromagnets. In this limiting situa- 
tion 7^ > — > while 7^ > is nonzero and independent 
of P. 

The quantum tunneling phenomenon leads to another 
interesting feature of the system under consideration. 
Generally, variation of the S z component of the SMM's 
spin by one is associated with spin reversal of one of 
the electrons in the leads. In the considered situation, 
an electron flips its spin from | to j, which is accom- 
panied by the transition of the molecule state from \m) 
to \m + 1). This process is continued until the probabil- 
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FIG. 3: (color online) The average charge Qlr transferred 
between the electrodes during the SMM's spin reversal pro- 
cess as a function of the field sweeping speed c for various 
polarizations of the right electrode P R . The polarization of 
the left electrode is P L — 0.4. The magnetic moments of the 
leads are kept either in the parallel (a) or in the antiparallel 
(b) configuration. The insets represent the DOS at the Fermi 
level in the leads, where the arrows indicate the direction of 
the charge transfer between the leads. 



ity of finding the SMM in the spin state \S) is equal to 
1, which corresponds to full reversal of the SMM's spin. 
Moreover, some of the spin-flip relaxation processes are 
associated with a transfer of a single electron charge from 
one electrode to the other. Owing to the spin asymmetry, 
the average charge transferred between the leads during 
magnetic moment switching becomes nonzero. As a con- 
sequence, one can expect a current pulse (I r ) associated 
with the spin relaxation after crossing each resonant mag- 
netic field [17] . The general formula for the relaxation 
current flowing from the left to right side of the junction 
is Ir = e££i El =1 Ps-m^T-l - 7^>), where e is 
the electron charge, and Ps- m as well as de- 
pend on the field H . Here, M is determined from the 
condition Ps-(m+i) ~ 0- The time scale on which re- 
laxation currents vanish is of the order of the relaxation 
times. For these reasons, current impulses are very short 
for the parameters corresponding to fig. [5] and would not 
be resolved. 

We note that the transition times only slightly change 
when the field is increased, and the corresponding time 
scale is considerably smaller than the time scale set 



by c. Therefore one can assume that they are con- 
stant, which allows us to find analytic expression for 
the total electronic charge transferred from the left to 
the right lead during the reversal of the SMM's spin, 
Qlr — e ^J2jLijPs-j, see figs. [3] T is a coefficient 
whose value depends on the magnetic configuration of 
the leads, and r j4P( p ) = 2(P L ± P R )/[A - {P L =f P r ) 2 ], 
respectively for the parallel (lower sign) and antiparallel 
(upper sign) magnetic configurations. 

The sign of Qlr, which is related to the direction of 
average charge flow, is determined by the magnetic con- 
figuration of the system. In the antiparallel configuration 
we observe only one direction of the average charge flow, 
i.e. from left to right. In the parallel configuration, on 
the other hand, the average flow can be either from left 
to right (P L > P R ) or from right to left (P R > P L ). For 
pL _ pR th ere i s n0 resultant flow of charge. The change 
of the charge transfer direction can be explained by con- 
sidering the DOS at the Fermi level in both electrodes 
and taking into account that a change in the SMM's spin 
state is possible only when it is associated with the flip 
of an electron from j" to J.. 

In conclusions, we have considered spin reversal of a 
SMM due to the phenomenon of QTM in a time depen- 
dent magnetic field and in the presence of two external 
spin-polarized leads. The molecule was assumed to be 
coupled to electrons in both leads via exchange coupling. 
We showed that full spin reversal is associated with a 
transfer of certain charge from one lead to the other, de- 
spite no bias voltage is applied. When external bias is 
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applied, then new effects may be observed [IE 
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